Physical aging and relaxation of residual stresses in a colloidal glass following flow 
cessation 
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Dilute Laponite suspensions in water at low salt concentration form repulsive col- 
loidal glasses which display physical aging. This phenomenon is still not completely 
understood and in particular, little is known about the connection between the flow 
history, as a determinant of the initial state of the system, and the subsequent aging 
dynamics. Using a stress controlled rheometer, we perform stress jump experiments 
to observe the elastic component of the flow stress that remains on cessation of flow 
or flow quenching. We investigate the connection between the dynamics of these 
residual stresses and the rate of physical aging upon quenching from different points 
on the steady state flow curve. Quenching from high rates produces a fluid state, 
G" > G' , with small, fast relaxing residual stresses and rapid, sigmoidal aging of 
the complex modulus. Conversely, quenching from lower shear rates produces in- 
creasingly jammed states featuring slowly relaxing stresses and a slow increase of the 
complex modulus with system age. Flow cessation from a flxed shear rate with vary- 
ing quench durations shows that slower quenches produce smaller residual stresses at 
short times which relax at long times by smaller extents, by comparison with faster 
quenches. These smaller stresses are correlated with a higher modulus but slower 
physical aging of the system. The characteristic time for the residual stress relax- 
ation scales inversely with the quench rate. This implies a frustrated approach to 
any ideal stress-free state that succinctly reflects the frustrated nature of these glassy 
systems. 
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I. INTRODUCTION 



Many complex fluids show some degree of time dependence in their properties. Such dis- 
plays of thixotropy are well known in c olloidal d ispers io ns for example , b oth under attractive 



and rep ulsive inter-particle potentials 



MewisI f ll979h : Barnesl fll997h : 



DuUaert and Mewis 



(j2005bl )]. When brought to rest, the properties of these systems m ay continue to evolv e 



as thermally dr i ven p ar ticle motion drives st ructural rear r angem ents 



Cipelletti et al. 



mm -. 



Abou et al. 



femil); 



Suarez et al 



Willenbacheil fllQQGl) : 



mm -. iNeri and Osuiil f l2009bh ]. 



a process known as physical aging. These rearrangements lead to slowly increasing viscosity 
and complex modulus as the system gradually "relaxes" and approaches some long time 
equilibrium or steady state at rest. There is qualitative similarity between the aging behav- 
ior of these colloidal systems and the glassy dynamics observed in polymer and molecular 
structural glasses , key among t hem th e display of time-elap s ed tim e rescaling of t he dy - 



namic properties 



Cloitre et al. 



fcoOOh : ICipelletti and RamosI (12002[ l: 



Cloitre et al. 



(12003[ l: 



Joshi and Reddyl ( 120081 )]. However, non-trivial distinctions still separate the behavior of 



McKenna et al. 



colloid al systems from that of traditional glasses, as highlighted recently 
( 120091 )]. The viscosity of these fluids under flow is comprised of both hydrodynamic and elas- 
tic or thermodynamic contributions. In dilute dispersions of Brownian hard spheres, based 
on Einstein's visc osity relation, the elastic and viscous terms have been given by Batchelor 



Batchelorl (119771 )] in Equations HE 



7/," = 0.97(^' 



(1) 
(2) 



where t]^ is the viscous contribution to the relative viscosity, 77^ is the elastic (Brownian) 

contribution and o? is the volume fraction of particles in the fluid. More recent models 
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have extended these considerations int o concentrated Brow nian dispersions |Bradyl (119931 



as well as weakly aggregated systems (Potanin et al. 



(I1995I )]. In flocculating systems, the 



elastic term depends on the strength of inter-particle attractions due to the presence of 



Sonntag and Russell (11987a 



(1987a, 


b); 


Osuji et al. 


(2008)] 



Deconvolution of the viscosity into rate dependent elastic and viscous terms can be made 
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via stress jump measurements where a suspension in steady state flow is subjected to a 
sudden drop to zero shear rate, 7 = 0. The viscous component of the shear stress dissi- 
pates nearly instantaneously, leaving only the elastic com ponent as the me a surab le quantity 



that contributes to the viscosity decay function, r]_ {t, 7) 



Watanabe et al. 



fll996h ]. Suitable 



modiflcation of strain controlled rheometers employing a force rebalancing transducer for 



torque measurements have enabled observa tions of the e^ 



as early as 20 ms a: 



ter cessation of flow 



Mackay et al. 



a stic s tre ss component a t time s 



fll992h : 



Kaffashi et al. 



fll997[ ): 



DuUaert and MewisI ( l2005al )] . Extrapolation of the time dependence of the stress back 



to the cessation of shear at t = allows the steady state flow curve to be represented in 
terms of the rate dependent elastic and viscous contributions to the system viscosity as the 
shear stress before the jump is given as o"(7) = a"" + a'^ = {1]^ + ri'^)'j. The difficulty of 
performing these measurements is matched by the insight they provide into the dispersion 
rheology. For example, in flocculating materials, the absence of hydrocluster induced shear 
thickening has now been attributed to the precipitous decrease in the elastic component of 
the shear stress with increasin g shear rate, compared to the r ather more modest increase of 
the hydrodynamic component Gopalakrishnan and Zukoskil (120041 )] . 

In addition to decoupling the flow stress components, stress jump measurements also 
permit the observation of the relaxation of the elastic stress. A single exponential form 
was observed in weakly floc culated colloidal dispersions of fumed silica and carbon black 



in hard sphere systems 



DuUaert and MewisI (l2005al)] out to ^ 0.1 s, whereas a power-law depend ence was found 



Kaffashi et al. 



(119971 1: lO'Brien and Mackavl (120001 )] out to ^ 1 s. 



Stress jump measurements have also been applied in polymer systems to follow the relaxation 

oxane) 



of deformed interfaces in immiscible blends of poly(isobutylene) in polyfdimethyls i 
(PDMS), where an exponential dependence was found, out to 4 s [Vinckier et all ( 1l997l )]. 
In more recent work on polymer-clay nanocomposites, single exponential relaxations of the 
elasti c flow stress component were observed, both for a Newtonian suspending polymeric 



f luid iMobuchon et all (120071 )] as well as for a low molar mass PDMS medium [Zhu et al. 



( l2009l )]. both out to ~ 0.2 s. Clearly, the dynamics of the elastic stress displayed on cessation 
of shear, the residual stress, so to speak, can provide important insights into the microscop- 
ies of the system under study. In materials which display physical aging, the same thermally 
driven structural rearrangements which give rise to increasing moduli are also responsible 
for the decay of the residual stress. Thus, one should expect that the time evolution of 



the system properties should be related in some way to the dynamics of this stress relax- 
ation. Surprisingly, this connection remains quite unexplored and residual stress dynamics 
in colloidal systems have typically been considered only at short times scales, versus the 
comparatively long times at which physical aging is studied. 

Here, we perform stress jump measurements on a Laponite dispersion and monitor the 
relaxation of the residual stress out to 100 s, as a function of the rate of flow arrest. We also 
study the physical aging of the system via small amplitude strain controlled dynamic tests 
following cessation of shear. Our data reveal striking parallels between the physical aging 
and residual stress relaxation dynamics. In particular, the characteristic time for the stress 
relaxation increases with the duration of the flow quench step. Concurrently, the rate of 
physical aging also decreases, implying a frustrated approach to the long time equihbrium 
state of the system. 



II. EXPERIMENT 
A. Materials 



Our system is an aqueous dispersion of Laponite which consists of disc-like clay particles 
of roughly 25 nm diameter and 1 nm thickness. Repulsive glasses are formed at low ionic 
strength in fairly di l ute co nditions due to long range electrostatic repulsion between the 



particles 



Bonn et al. 



Knaebel et al. 



cal means 



mm-. 



1998(). The svs t em di spl ays glassy dynamics, as studie d by scattering 



Bonn et al 



Abou et al. 



(12002); 



Abou et al 



(120011) : iBandvopadhvav et al\ (2004f)] and rheologi- 



Joshi and Reddy 



mm ]. Samples of 



((20031); 

concentrations = 3, 4 and 5 wt.% are prepared by dispersion of Laponite XLG powder 
(Southern Clay Products) into ultra-pure w ater, adjusted to pH 9 .5 via addition of NaOH 



to ensure chemical stability of the particles 



Mourchid and Levitzl ( ll998l )]. The samples are 



agitated vigorously for 20 minutes and then allowed to develop completely for 5 days under 
quiescent conditions, defining a well controlled and reproducible initial state. 



B. Method 



Rheology is conducted in strain controlled mode using an MCR301 instrument (Anton- 
Paar) in the cone-plate geometry (1°, 50 mm diameter steel cone). Evaporation of water 
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from the sample was successfully suppressed without perturbing our system via application 
of a thin film of mineral oil at the sample edge. Steady state flow curves were collected from 
7 = 10^^ to 10^ s^^ using a logarithmically varying equilibration time from 60 s at the lowest 
shear rates to 1 s at the highest shear rates. Conventional stress relaxation data were taken 
after a resting period of 300 s following a pre-shear at 100 s^^, using a strain of 7 = 3%, 
which is within the linear viscoelastic regime for all compositions studied. Frequency sweeps 
were carried out using this same strain from 100 to 1 rad/s. 



1. Variation of pre-shear rate 

Samples were subjected to varying pre-shear flows at rates between 4 and 1000 for 
a duration of 100 s. In each case, this time was more than sufficient to establish a robust 
steady-state in which the viscosity of the sample was well defined within ±5%. The samples 
were subjected to a linear quench of 1 s duration (tg = 1 s) from the fiowing state to the 
stationary or at-rest state. Physical aging and stress relaxation measurements were carried 
out on the system following separate pre-shear and mechanical quenching steps as just 
described. The physical aging of the samples after cessation of fiow was monitored via the 
time evolution of the complex shear modulus in the linear viscoelastic regime using strain 
controlled oscillations at 7 = 3 % at an angular frequency of = 10 rad/s. The relaxation 
of the residual elastic component of the fiow stress in the system after cessation of fiow, 
lere termed the residual stress, was then characterized, as described in a prior publication 



Negi and Osujil ( l2009cl )]. Briefiy, a constant zero strain rate condition was imposed on the 
sample. The stress required to maintain this 7 = stationary condition is the residual 
stress. Checks were performed using a Newtonian fiuid, which should have no measurable 
residual stress, to verify the ability of the instrument to accurately make the measurements 
reported here. 



2. Variation of quench duration 

Samples were pre-sheared at a fixed rate of 7 = 100 s~^ for 100 s and then brought to 
rest via a linear ramp over varying durations, tg, from 1 to 300 s. Following this cessation of 
fiow, physical aging was monitored as above via the time evolution of the complex modulus 
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in the linear viscoelastic regime, 7 = 3.0% and u = 10 rad/s. Following the same pre-shear 
and flow cessation, the relaxation of residual stresses were monitored, as described above. 



III. RESULTS AND DISCUSSION 
A. Effects of Pre-Shear Rate 



Laponite forms a yield stress fluid when dispersed in water, even at relatively low con- 
centrations. As a result, the steady flow curves do not display Newtonian plateaus at low 
shear rates, as shown in Figure [TJ The fluids display monotonic shear thinning, with shear 
thinning exponents |?t,| > 0.9 for all three concentrations considered here. We do not ob- 
serve high shear rate viscosity plateaus, even out to 10^ s^^. This strong she ar thinning 



behay ior over a wide range of shear rates is typical of Laponite suspensions 



Bonn et al. 



( 120021 ): 



Abou et al 



( I2OO3I )]. Additionally, we do not observe any stress plateaus in the 
flow curves, suggesting that there are no obvious intervening regimes of shear banding or 
other flow-induced structuring. Absence of such makes these suspensions ideal candidates 
for studies as conducted here. 




0.01 0.1 



1 10 100 

shear rate (1/s) 



1000 10000 



FIG. 1. Flow curves for 3, 4 and 5 wt.% samples on decreasing shear rate. 



Stress relaxation measurements on cessation of flow were performed for shear rates be- 
tween 4 and 1000 s~^, as shown in Figure El The data show a slowing rate of decay of 
the residual stress as the initial shear rate of the system is decreased. In some cases, at 
longer times, there are slight increases in the stress which are attributed to structuring of 
the system at rest, which causes the viscosity of the suspension to rise with time. Because 
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the rheometer in practice cannot maintain an identically zero shear rate (in general the 
minimum controllable was 7 ~ 10^^ s^^), this results in an increasing component of the 
stress that may become measurable at long times. The data for 3 wt.% following cessation 
from 7 = 300 s^^ are particularly egregious however, for reasons that remain unclear. The 
instantaneous value of the residual stress immediately upon flow cessation is required to 
accurately determine the elastic component of the flowing stress at different points on the 
flow curve. Our measurements, however, did not permit observation of the residual stress 
at times shorter than 0.1 s. We take ar{t = 0.2 s) as representative of the elastic stress 
component, o"^, during the immediately preceding flow, rather than performing a log-linear 
regression of the data to t = 0, which, strictly speaking is required for this purpose. Per- 
forming a log-linear regression would imply the assumption that the same functional form is 
maintained at shorter times, which is something that we cannot justify with complete cer- 
tainty. Despite the quantitative difference between the two estimates, the data do provide 
the same qualitative result. The elastic stress component can be expressed as a fraction 
of the overall flow stress as R'^i'j) = crr{i',t ~ 0)/'^(7)- As shown in Figure [3l the elastic 
stress component decreases dramatically with increasing shear rate, illustrating the strong 
reduction in the cohesion of particle assemblies under flow. While this approach results in 
an underestimate of the elastic stress component, it does capture the large decrease in 
with increasing shear rate that is expected in these systems. 

Physical aging measurements following rate quenches from the same flow rates as above 
show that at long times, all samples followed what appeared to be a weak power-law evolution 
of the complex modulus with time. Figure HI However, the behavior at short times varied 
strongly with the initial shear rate, and the vertical separation at long times increased with 
increasing volume fraction. Quenches from the lowest shear rate, 7 = 4 s~^ produced a higher 
complex modulus, with G' > G", at short times which evolved with a monotonic near power 
law display throughout. By contrast, quenches from the highest shear rates 7 = 1000 
produced systems that were initially fluid-like, with G" > G' at short times, for = 3 wt.%. 
The complex modulus evolved in a sigmoidal fashion, rapidly increasing at intermediate 
times before assuming the weak near power law growth at long times. Quenches from 
intermediate shear rates showed a smooth variation between these two limiting behaviors. 
No intersection of the different trajectories was observed during the measurement. 

The sigmoidal variation of the modulus in time is quite striking, but not entirely unex- 
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pected. While many dis ordered colloidal sys te ms display power-law ag inK as implicit in time 



elapse d time rescaling [Cloitre et all ( l2000l ) 



Knaebel et al. 



mm -. iRamos and Cipelletti 



( 120051 )]. intuitively, one understand that the power law dependence cannot persist indefi- 
nitely, if it is assumed that the system actually has a well defined, finite state or modulus 
at long time. Part of the difficulty in studying these systems, however, is the slow nature of 



8 



I I I I I I 



10 100 1000 



Preshear Rate (1/s) 

FIG. 3. Ratio of elastic stress at t=0.2 s to the overall flow stress as a function of the shear rate. 
Lines are drawn as a guide to the eye. 

the dynamics, and the slowing of the dynamics with the approach to equilibrium or steady 
state. Here we see that by starting from a highly fluid state characterized by a vanishing 
elastic stress component, we are able to access the short time regime of the system where 
it displays a sigmoidal response as it ages at first slowly, then rapidly, then slowly again as 
it approaches its long time steady state. This rapid evolution of the complex modulus is 
mirrored in the rapid relaxation of the residual stress for samples quenched from the highest 
shear rates. For a fixed shear rate, moving to higher concentrations increases the residual 
stress component and slows the dynamics of the system overall. A Carreau-like model, 
Equation [3] was used to estimate a characteristic relaxation time for the sigmoidal evolution 
of the complex modulus. 



G*{t) = G*{oo) + 



G*(0) -G*(oo) 



(3) 



[1 + it/ry] 

The rate of aging can be represented by the parameter p which gives an indication of the 
steepness of the curve near the midpoint, t = r, of the growth between presumed constant 
t = and t = oo values. The inset plots of Figure H] show the increase of p with increasing 
shear rate before flow cessation, as described. Interestingly, a sigmoidal response has also 
been observed in the time evolution of the elastic modu lus of a system composed of multi-arm 
star polymers in solution |Christopoulou et all (120091)] The soft colloidal glass thus formed 



has important differences from the currently considered Laponite sus pension. Notab 



in the polymer glass proceeds via a complex two-step mechanism 



Helgeson et al. 



V ; agin g 
J2OO7I ): 



10 



100 



1000 




t(s) 

FIG. 4. Evolution of the complex modulus for 3, 4 and 5 wt.% samples (top to bottom) following 
stress jumps from different shear rates on the flow curve as indicated. The lines are fits to the data 
using Equation [3l Oscillations were performed using a strain of 7 = 3 % at an angular frequency of 
10 rad/s. Inset: Value of the exponent p in Equation [3] describing the slope of G*{t) in the vicinity 
of t = r. Lines are drawn guide to the eye. 



Rogers et al\ ( 120101 )]. and the system is able to access a presumed equilibrium state at 
long times (~ 5E4 s) where the modulus is stable in time. Nevertheless, the parallels are 
noteworthy. 

The correlation between the earliest measured modulus of the system and the short time 

10 



value of the residual stress is shown in Figure |5l Smaller residual stresses are produced 
by quenches from higher shear rates and are correlated to lower values of the complex 
modulus. A linear dependence has been observed between the elasti c modulus a nd th e 



Osuji et al. 



mm -. 



residual stress in par t iculat e dispersions with attractive interactions 
Ovarlez and Chateaul (120081 )]. Here, the dependence is roughly linear, although for reasons 
that remain unclear, the low stress data for (y9 = 4 wt.% are somewhat non-conforming. 
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FIG. 5. Relation between the modulus and the stress on shear cessation, for different pre-shear 
rates. 



B. Effects of Quench Duration 

The preceding data clearly show that the dynamics of physical aging are paralleled by 
those of residual stress relaxation. This is not entirely surprising as the same thermally 
driven particle motion that gives rise to rearrangements that age the system also result in 
the relaxation of the elastic stresses built into the system on cessation of flow. Additionally, 
aging from a more fluid state should proceed faster than from a more jammed state. If 
we now consider quenches from a fixed shear rate, but over varying durations, we should 
expect to observe at some point the finite effects of a non-steady state flow history on the 
sample. A rapid flow cessation transitions the system from a free flowing steady state to a 
stationary, jammed state in a short time. Increasing the duration of the quench allows the 
system to sample configurations at lower shear rates en-route to the stationary state. The 
equilibrium condition at any particular shear rate is provided by the known viscosity from 
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the steady state flow curve and so the instantaneous viscosity during the quench provides 
a good metric for the departure of the system from equihbrium during the cessation of 
flow. Over a range of quench durations, the system should display a spectrum of behavior 
in which the viscosity at any given shear rate during the quench deviates from the steady 
state value by an amount that scales with the quench rate. We monitor the instantaneous 
viscosity during the cessation of flow and find precisely this behavior, shown in Figure |6l 
As expected, the magnitude of the shear thinning exponent, n where 77(7) ~ 7" decreases 
continuously, in a linear fashion, as the quench duration is decreased. 

100000 

> 

"5) 

8 10000 

at 
> 

> 
« 

1000 

1 10 100 

shear rate (1/s) 

FIG. 6. Different quench rates result in the samples traversing different paths in the ergodic to 
non-ergodic transition on cessation of flow. Data are for a 4 wt.% sample. Inset: Shear thinning 
exponent during the quench decreases in magnitude as a function of decreasing quench duration. 

Prior to measurement of residual stresses, as well as for the modulus measurements, 
samples were subjected to a rejuvenating flow at the selected shear rate, 7 = 100 s~^. 
During these rejuvenating steps, the flow stress and thus the viscosity evolve to well defined 
levels that are characteristic of the steady state of the system at the applied flow rate (insets. 
Figure [7]). The elastic stress that results on flow cessation is due to the structuring of the 
system that takes place during the finite quench interval. Surprisingly, we find that fast 
quenches lead to a higher residual stress than produced for slow quenches. Additionally, 
the residual stresses relax faster and to lower values for fast quenches than those produced 
by slower quenches. Figure [71 The data are strikingly different in form than those for a 
conventional stress relaxation, in response to a step-strain. They exhibit neither a power 
law nor stretched exponential dependence on time, but instead are suggestive of a sigmoidal 
approach to a finite, non-zero value at long times, much like the residual stress data of 
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Figure [2l Notably, the relaxation shows a clear and monotonic dependence on the quench 
duration, tq. 
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FIG. 7. Evolution of residual stress for different quench durations as indicated, for 3, 4 and 5 wt.% 
samples (top to bottom) quenched from 100 s~^. Lines are fits to the data using Equation HI The 
insets show the flow stress during the rejuvenation steps, indicating that the system starts from a 
consistent state. 

We fit our data with good fidelity again using a sigmoidal function. Equation H] for the 
time dependence of the residual stress, ar{t). Here, again, r is a characteristic time that 
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FIG. 8. Quench time dependence of the characteristic stress relaxation time, r derived from 
Equation m The hne shows a slope of 1. 

describes the decrease of the sigmoidal function to its vertical midpoint between the value at 
t = 0, (7(0) and at t = oo, cr(oo), and the exponent p characterizes the slope of the function 
in the vicinity of t = t. Remarkably, we find that the characteristic time for the relaxation 
of the residual stress is a linear function of the quench time, r ~ tg, as shown in Figure [SI 



The modulus of the glass, measured immediately after cessation of shear is inversely 
proportional to the rate of quench. The slowest quench rates, with flow cessation over 
periods of 300 s, lead to the highest moduli. The modulus shows a slow increase that is well 
described by a power law at late times, but with significant sigmoidal character at shorter 
times. Figure M This is not unlike the response observed for quenches from different initial 
shear rates. Attempts to superpose the data by time-elapsed time rescaling fail. Figure 
[TOl That is, shifting in time to account for the true sample age due to passage of time 
during the quench does not produce a satisfactory superposition of data. These results 
show that the flow history of the sample is important in determining its aging path, and 
different trajectories through the non-steady state flow curves en route to the stationary 
state invariably result in different and non-superposable aging behaviors thereafter. The 
correlation between modulus and residual stress is shown in Figure [TH Here now, there is 
no longer a simple direct correspondence between stresses and moduli. The residual stress 
displayed is not just function of structuring that occurs during the flnite quench, but also 




(4) 
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relaxation of elastic components that were supported at the start of the quench. 
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FIG. 9. Evolution of the complex modulus in 3, 4 and 5 wt.% samples (top to bottom) quenched 
from 7 = 100 s^^. Oscillations are performed using a strain of 7 = 3 % at an angular frequency of 
10 rad/s. The legend indicates the duration of the quench, tq, in seconds. 

The data for the residual stress relaxation can be contrasted with that obtained by tradi- 
tional step-strain experiments as well as frequency dependent moduli. Frequency sweeps of 
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FIG. 10. Evolution of the complex modulus, shifted according to the quench durations. 
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FIG. 11. Relation between the modulus and the stress on shear cessation, for different quench 
rates. 

the suspensions show a very weak dependence of both storage and loss moduli on frequency, 
Figure [121 Data was collected after quenching from 7 = 100 s~^ over 1, 10 and 100 s, follow- 
ing which the samples were allowed to sit quiescently for 300 s prior to measurement. This 
ensured that the measurement time (~ 100 s) was small enough compared to the sample age 
that changes in the system properties due to aging during the measurement were reasonably 
small. Similarly, step strain experiments were performed after an 300 s quiescent period to 
evaluate the conventional stress relaxation behavior. A strain of 7 = 3% was used, which 
is in the linear regime for the system, as confirmed by strain sweep experiments. The data. 
Figure [131 are in good agreement with the slow relaxation that is inferred from the G{u) 
data of Figure [121 Less than half a decade of relaxation is observed over 3 decades of time 
for 3 wt.%, and substantially less than this for 4 and 5 wt.%. For both G{uj) and G{t), there 
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FIG. 12. Frequency sweep measurements. Samples were quenched over the 3 different times 
indicated and aged for 300 s before measurement. 

is little dependence on the quench rate both in the absolute values of the moduli measured, 
and, more significantly, in their time(frequency) dependences. 
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FIG. 13. Conventional stress relaxation measurements. Samples were quenched over the 3 different 
times indicated and aged for 300 s before measurement. 




IV. CONCLUSION 

In summary, our investigation has shown that the quench rate, as well as the initial 
flow rate from which the system is quenched, both strongly influence the aging dynamics. 
The effects of the initial flow rate are particularly intuitive, assuming an inflnitely fast 
quench such that the properties (structure) immediately on cessation of flow are those of 
the immediately preceding flow regime. One can readily imagine that a system that is in 
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steady flow at an infinitely low shear rate, 7 — )■ 0, would display an infinitesimal rate of 
aging upon cessation of flow. Conversely, a system in steady flow at a very high shear rate, 
7 — )■ 00, is completely fluidized, and should display the the full aging response from fluid 
to semi-solid as the system undergoes structural arrest and further vitrifies with passage of 
time after cessation of flow. For all finite initial shear rates, the displayed behavior should 
fall on a spectrum between these two limiting described in Figure [HI Whether or 

not a power-law or logarithmic or even exponential aging response is observed then becomes 
a function of the initial flow rate of the system and the duration of the observation window, 
among other fa ctors. Promising l y, rece nt work in this a rea has been su c cessfu l in revealing 



both the stress 



Negi and Osujil ( l2009al )] and frequency jNegi and Osujil ( I2OIOI )] dependence 



of the complex modulus during structural arrest following rapid cessation of flow from high 
shear rates. The behavior observed for the currently studied Laponite system is in broad 
agreement with this framework of a globally sigmoidal response which we advance as a 
generic response for out of equilibrium materials. This argument has also been given based 
on earlier work in polymer glasses, where a sigmoidal form is suggested for a perfect thermal 



quench, with power - 



window 



aw a ging simply the consequence of a limited observation or aging time 



McKennal ( 120031 )]. Here, abrupt mechanical quenching of the flowing suspension 



permits access to a larger spectrum of the system response. The lack of superposition 
of the complex modulus with waiting time may reflect the different energy landscapes or 
distribution of relaxation mechanisms available to the system as a function of the initial flow 
rate and the fact that the evolution of the landscapes is not along a single trajectory, but 
depends on the initial condition from which the system is evolved. 




FIG. 14. Generic picture of the evolution of the complex modulus. 
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With a fixed initial condition, the quench rate is seen to exert a strong influence on the 
aging dynamics in this repulsive colloidal glass and serves to characterize the departure of 
the system from equilibrium. Rapid quenches lead to faster aging from smaller moduli. 
Further, the aging dynamics are well correlated to the decay of residual stresses that are 
trapped in the glass at the ergodic to non-ergodic transition on cessation of flow. These 
stress dynamics are dependent on the quench rate and display a linear dependence of their 
characteristic relaxation time on the duration of the quench. Notably, they follow neither 
a simple exponential nor a power law form. The simultaneously decreasing initial value 
and final extent of relaxation of the residual stress concisely represent the frustrated and 
continuously slowing approach to equilibrium in these materials. This work suggests that 
the flow cessation rate may be identified as somewhat akin to an effective cooling rate 
for the colloidal glass wherein infinitely slow cooling, — )■ oo, maintains equilibrium and 
results in a stationary state on complete flow cessation. Many of the features displayed here 
resonate intuitively, and should be observed in other colloidal systems which display glassy 
or glass-like dynamics. Dynamic scattering can provide useful data in these systems and 
measurements which incorporate control over the flow arrest are a natural follow on to this 
work. 
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